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Abstract 

An exceptionally stable hollow cage containing 20 scandiums and 60 carbons was iden¬ 
tified. This SC 20 C 60 molecular cluster has a 7), point group symmetry and a volleyball-like 
shape, that we refer to below as "Volleyballene". Electronic structure analysis shows that the 
formation of delocalized 7t bonds between Sc atoms and neighboring five-membered carbon 
rings is crucial for stabilizing the cage structure. A relatively large HOMO-LUMO gap was 
found. The results of vibrational frequency analysis and molecular dynamics simulations also 
demonstrate that this Volleyballene molecule is exceptionally stable. 

Since the experimental observation of C(,q,— many very interesting structures have been pro¬ 
posed, such as M$C \2 (M=Ti, V, Zr, Hf, et al)2£ Au 2 o,- AU 32 ,- AU 42 ,- M@Si w (M=Transition 
Metals; «=14,15,16),- Eu@Si 2 o,- Eu 2 Si 3 o,- Bgo,— and B 40 .— In the present work, an excep¬ 
tionally stable hollow cage, SC 20 C 60 Volleyballene, has been identified by first-principles density 
functional theory studies. 

*To whom correspondence should be addressed 
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Our first-principles calculations were performed within the framework of spin-polarized den¬ 
sity functional theory (DFT). The exchange-correlation interaction was treated within the general¬ 
ized gradient approximation (GGA) using the Perdue-Burke-Ernzerhof (PBE) exchange-correlation 
functionals. 12 A double-numerical polarized (DNP) basis set 13 was chosen to carry out the elec¬ 
tronic structure calculation with unrestricted symmetry. For the transition metal atoms, relativistic 
effects in the core were included by using the DFT semi-core pseudopotentials (DSPP).— 
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Figure 1: The configurations of the Volleyballene SC 20 C 60 viewed in (a) ball and stick model 
(Large ball: Sc atom; small ball: C atom) and ( b ) CPK style. Part (c) shows a volleyball. 

Figure 1 shows the configuration of the Volleyballene Sc2oC6o- It is comprised of 20 scandium 
atoms and 60 carbons with a volleyball-like shape. This new, lowest energy SC 20 C 60 molecule had 
T h point group symmetry and robust stability. As may be seen from Fig. 1, it may be viewed as 
consisting of six ScgQo subunits joined together in a crisscross pattern. In this structure, there are 
12 pentagonal rings made of carbon atoms (C-pentagons) and 6 octagonal rings of scandium atoms 
(Sc-octagons). Every group of two C-pentagons is surrounded by one Sc-octagon to give a Sc^Cio 
subunit. 

The 20 Sc atoms link to form 12 suture lines, and the average Sc-Sc length is 3.222 A. De¬ 
pending on the coordination of the Sc atoms, they may be divided into two types, Sc 1 and Sc 11 (see 
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Fig. 2b). For the C-pentagons, there are three C-C double bonds and two C-C single bonds. Of the 
three C-C double bonds, one has a length of 1 .466 A and the other two both have lengths of 1 .434 

o o o 

A. The two C-C single bonds are both 1.446 A long. Along with a 1.463 A C-C bond connecting 
the two C-pentagons, the average C-C bond length is found to be 1.446 A. The average Sc-C bond 
length is 2.248 A. 

The following aspects of the stability of the Volleyballene SC 20 C 60 were investigated: the bond¬ 
ing character, the binding energy, the vibrational frequencies and the molecular dynamics, the latter 
through ensemble simulations. 

Firstly, the bonding character of Volleyballene SC 20 C 60 was investigated by analyzing its de¬ 
formation electron density, as shown in Figs. 2b and 2c. For the C atoms, there are obvious 
characteristics of .syr-like hybridization, and each C atom has three <7 bonds. For the C atoms 
that neighbor pairs of Sc atoms, two of the three orbital lobes point towards the two neighboring 
C atoms, and the third lobe points towards the center point of the line joining the two Sc atoms 
but with a slight inclination toward Sc 1 . For the other cases, the C atoms have three neighboring 
C atoms and form a bonds between all four C atoms. An outer delocalized electron is then left 
over for each carbon atom. The Sc atom has three outer electrons, one lying in the 3d state and 
the other two lying in the 4.s state. In each So^Qo subunit, a delocalized n bond forms between 
the two Sc 1 atoms and one C atom located at the waist position. This is crucial for stabilizing the 
Scj^Cio subunit. For the Sc 11 atoms, there is a chiral n bond three petals pointing towards the three 
neighboring C-C bonds, which strengthens the link between the Sc^Cio subunits. 

In order to further validate the stability of the Volleyballene SC 20 C 60 , we constructed five other 
molecular structures with different combinations of 12 C-pentagons and 20 Sc atoms. Calculations 
were carried out within the same framework as described above. After energy minimization, it was 
found that the optimized structures all showed relatively large deformations and did not retain their 
original topologies. However, closer observation of these configurations indicated that although the 
overall configurations were not viable, one or more ScgCio subunit usually appeared, which further 
suggests that the ScgCio subunit is very stable. Figure 3 shows three of the cage-like geometries 
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Figure 2: (a) PDOS for the Volleyballene SC 20 C 60 , (b,c) the deformation electron density viewed 
from the top of Sc^Cio subunit and the suture line (c), and (cl) HOMO, (e) LUMO orbitals for 
Volleyballene SC 20 C 60 • The zero of energy for the PDOS curve (a) is taken to be the Fermi 
energy, and the length of each vertical line represents the corresponding relative amplitude for 
each electronic orbital. The isosurface for the deformation electron density (parts b and c) and the 
two orbitals (parts d and e) are taken to be 0.030 and 0.005 e /A 3 , respectively. 
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examined before and after optimizations. Beneath each isomer is listed the relative binding energy 
(A Eb) with respect to the Volleyballene. Clearly, from the energy point of view, the Volleyballene 
(shown in Fig. 1) was the lowest energy structure. 



Figure 3: Three kinds of typical initial and optimized configurations of SC 20 C 60 clusters. Beneath 
each isomer is listed the relative binding energy (AE h ) with respect to Volleyballene SC 20 C 60 • Key: 
large ball, Sc atom; small ball, C atom. 

Next, a vibrational frequency analysis was carried out for the lowest energy Volleyballene, 
Sc2oC6o- There were no imaginary frequencies and the two highest intensity frequencies were 
found to be 468.9 and 472.3 cm -1 . Figure 4 gives the calculated Raman spectrum with a 300 K 
temperature and 488.0 nm incident light, in order to simulate a realistic Raman spectrum that can 
be compared to experimental results. The calculated results show that the Volleyballene SC 20 C 60 
is kinetically the most stable isomer. 

In addition, ab initio molecular dynamics simulations with NVE and NVT canonical ensembles 
were carried out to test the thermodynamic stability of the SC 20 C 60 Volleyballene. For the NVE 
ensemble, the total simulation time was set to be 2.0 ps with a time step of 1.0 fs at initial tem¬ 
peratures of 1000, 1400, 2000, and 2400 K. The results of the NVE simulations showed that the 
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Figure 4: The simulate Raman spectrum for the Volleyballene SC 20 C 60 with a 300 K temperature 
and 488.0 nm incident light. The Lorentzian smearing is set to be 20.00 cm -1 . The labels present 
the frequencies corresponding to the peaks of the intensities. 

structure of Volleyballene was not disrupted over the course of a 2.0 ps dynamical simulation at an 
initial temperature of 2400 K, equal to a ~1200 K effective temperature. For the NVT ensemble, 
the Gaussian thermostat— was chosen, and the total run time was set to be 1.0 ps with a time step 
of 1.0 fs at temperatures of T = 800 and 1000 K. The NVT dynamical simulations indicated that 
the Volleyballene SC 20 C 60 also retained its original topological structure up to a temperature of 
1000 K. These results indicate that the SC 20 C 60 Volleyballene has good thermodynamic stability. 

Considering further the question of stability, it is natural to explore the electronic structure. 
To this end, we calculated the partial density of states (PDOS) and frontier molecular orbitals, 
including the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular 
orbital (LUMO) of Volleyballene as shown in Figures 2a, 2d, and 2e. From the contours of the 
HOMO and LUMO orbitals, it may be seen that the HOMO orbitals are mostly localized on the 
C atoms and Sc 1 atoms. There is also obvious p-d hybridization. As for the case of the LUMO, 
the energy level features show that the LUMO orbital is doubly degenerate. On the Sc atom, there 
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are obvious d-orbital characteristics, with <7,2-like orbitals for Sc 11 atoms and other <7-like orbitals 
for Sc 1 . The LUMO orbital hybridization is predominantly sp-d hybridization. Close examination 
of the PDOS (see Fig. 2a) further confirms the hybridization characteristics of the HOMO and 
LUMO orbitals. All these results demonstrate that the hybridization between Sc d orbitals and C 
p orbitals is essential for stabilizing the cage structure. For the Volleyballene SC 20 C 60 , a relatively 
large HOMO-LUMO gap of 1.471 eV was observed, mainly due to the energy of the d atomic 
orbitals being much lower than that of the p orbitals. The Volleyballene SC 20 C 60 should therefore 
be an exotic buckyballene variant with exceedingly high chemical stability. 
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